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Abstract 

Multiple-Input Multiple-Output (MIMO) based Medium Access Control (MAC) protocols have 
received a good deal of attention as researchers look to enhance overall performance of Ad Hoc networks 
by leveraging multi antenna enabled nodes |[ll-|l5l. To date such MAC protocols have been evaluated 
through comparative simulation based studies that report on the number of concurrent links the protocol 
can support. However, a bound on the maximum number of concurrent links (MNCL) that a MIMO 
based MAC protocol should strive to achieve has hitherto been unavailable. In this paper we present a 
theoretical formulation for calculating the bound on the MNCL in a Mobile Ad Hoc Network (MANET) 
where the nodes have multiple antenna capability, while guaranteeing a minimum Quality of Service 
(QoS). In an attempt to make our findings as practical and realistic as possible, the study incorporates 
models for the following PHY layer and channel dependent elements: (a) path loss and fast fading effects, 
in order to accurately model adjacent link interference; (b) a Minimum Mean Squared Error (MMSE) 
based detector in the receiver which provides a balance between completely nuUing of neighboring 
interference and hardware complexity. In calculating the bound on the MNCL our work also delivers 
the optimal power control solution for the network as well as the optimal link selection. The results are 
readily applicable to MIMO systems using Receive Diversity, Space Time Block Coding (STBC), and 
Transmit Beamforming and show that with a 4 element antenna system, as much as 3 x improvement in 
the total number of concurrent links can be achieved relative to a SISO based network. The results also 
show diminishing improvement as the number of antennas is increased beyond 4, and the maximum 
allowable transmit power is increased beyond 10 dBm (for the simulated parameters). 
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I. Introduction 

Wireless Ad Hoc networking has emerged as an important aspect of next generation communi- 
cation systems. For conventional Single-Input Single-Output (SISO) system, interference among 
nodes drastically limits the number of concurrent (simultaneous) links in Ad Hoc networks. 
Multi antenna, multi-input-multi-output (MIMO), based wireless communications has the ability 
to spatially null interference and in so doing increase the number of concurrent links within a 
Mobile Ad Hoc Networks (MANET), thus increasing overall network throughput. In fact some 
work found in the literature [[Il-O look to MIMO capable MAC protocols as a means of 
increasing the network efficiency and its sum- throughput. 

The maximum number of concurrent links is a metric used in the literature [[T|, [|71, [HI 
to evaluate the capacity of a network. Examples of MACs that support concurrent links in 
a network where all nodes have multiple antennas at their disposal can be found in [[I]|-[l5]|, 
[|71. For convenience these MAC protocols will be referred to as concurrent-based MACs in 
this paper. The Null-Hoc [|2|| and SPACE-MAC [3] protocols look to enable concurrent links 
by using the Gram-Schmidt Orthonormalization, so as to create orthogonal channels among 
different links. The protocol in H uses adaptive interference cancellation both at the transmitter 
and at the receiver, as well as a simple power control protocol for each link. Multiple links 
are assumed to access the channel sequentially and work simultaneously. The MIMA-MAC 
protocol [l5| uses space division multiple access techniques to guarantee the concurrency of 
different communicating links in the network. 

Although these concurrent-based MAC protocols have proved to outperform the conventional 
SISO based MACs such as the IEEE 802.11 DCF [9], a natural question to ask is that how 
close they actually come to the theoretical bound (limit) of concurrency in Ad Hoc networks. 
Furthermore, since MIMO systems enable a variety of approaches in utilizing multiple antennas 
in the physical layer [IQ], it is also concerned that how this bound of concurrency is affected by 
the choice of MIMO algorithms and associated physical layer techniques. In this work we identify 
the theoretical Maximum Number of Concurrent Links (MNCL) in the network by considering 
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the following PHY layer and channel dependent elements: (a) path loss and fast fading effects; 
(b) different MIMO transmit/receive algorithms; (c) a Minimum Mean Squared Error (MMSE) 
based detector in the receiver; (d) optimal power control and optimal link selection. The derived 
MNCL acts as a performance benchmark for concurrent-based MAC protocols, and is also used 
as a metric for comparing different MIMO techniques and parameters. 

Our study is based on the assumption that each transmit/receive pair requires the same Constant 
Bit Rate (CBR). In this way, the MNCL that can be had in a MIMO capable network is identified 
subject to a minimum Quality of Service (QoS) constraint. In our case the QoS constraint is 
the received Signal to Interference plus Noise Ratio (SINR) which is directly related to Bit 
Error Rate (BER) and Packet Error Rate (PER). The proposed framework is executed via an 
iterative process for power allocation and a Backtracking-based search strategy for link selection. 
Results of different MIMO transmit algorithms such as Receive Diversity, Space-Time Block 
Coding (STBC) and Transmit Beamforming on the MNCL bound are studied. Using this bound 
we present design-relevant insight regarding the impact of the number of pairs, the number of 
antennas, and the maximum allowable transmit power per pair on the MNCL. 

The reminder of this paper is organized as follows. The proposed system model is described in 
Section II. The definition of concurrent links is provided in Section III. Optimal power allocation 
is investigated in Section IV. In Section V, we propose a backtracking-based strategy to find the 
optimal link selection. Numerical results are shown in Section VI, and we conclude in Section 

vn. 

II. System Model 

A. Network Model 

We focus on a small network (or sub-network), where each node is within the transmission 
range of any other node subject to path loss and Rayleigh flat fading (Fig. 1). Assume K 
simultaneously communicating pairs in the network. Each pair is composed of one transmit node 
and one receive node, which are all randomly distributed in the network. Nodes in different 
pairs are unique and independent. They are all equipped with M receive antennas, sharing 
the same frequency band and requiring the same Constant Bit Rate (CBR). The number of 
transmit antennas varies depending on the MIMO technique being considered. For simplicity. 
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each communicating pair is labeled as a Transceiver Pair. For the A;th transceiver pair, the 
transmit node and receive node are named Tx Node k and Rx Node k, respectively. 

The propagation between nodes is characterized by path loss and Rayleigh fading [fTT]| - [fT3ll . 
For path loss, we use the simplified model in |fT3l , which is: 

Lp{d){dB) = Lp{do) + lOalogio j- (1) 

"0 

In this model, a is the path loss exponent, do is the reference distance, and d is the distance 
between nodes obtained from node topology information. Both a and d^ are parameters in our 
study and can be set by the user. We use pkj to denote the power loss ratio from Tx Node j to 

Lp(d) 

Rx Node k, and p^j = 10 lo . 

We assume a flat fading environment, which is modeled by a Rayleigh distribution. Hfcj(m) 
denotes the Rayleigh fading channel from the mth antenna of Tx Node j to Rx Node k. It is 
an M X 1 vector and consists of independent identically distributed (i.i.d.) complex Gaussian 
random variables. 

The background white noise is a circularly complex Gaussian vector with covariance matrix 
(J%Im- is an M X M unitary matrix, and cr^ is given by: 

a^(dBm) = Vn + lO logio ( VT) + (2) 

W^(Hz) is the bandwidth of the system, while r7„(dBm/Hz) and F„(dB) are the power spectral 
density of the thermal noise and the noise figure of the receiver, respectively. They are assumed 
to be identical for all the nodes. 

B. Transceiver Model 

Three different MIMO techniques are considered in this paper, they are: (a) 1 x M Receive 
Diversity; (b) 2 x M Space-Time Block Coding (STBC); and (c) M x M Transmit Beamforming. 
In this paper, we assume that data symbols in the baseband have the same modulation type, 
regardless of single or multiple antennas. For fair comparison, the number of transmit antennas 
in these MIMO techniques are selected so that their spectral efficiency is the same as that of a 
Single-Input Single-Output (SISO) system. 
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1) 1 X M Receive Diversity: In the 1 x M scenario [fT4l . only one transmit antenna is employed. 
Without loss of generality, the 1st transmit antenna is used. Considering the k\h transceiver pair, 
the received signal at Rx Node k is given by: 

Yfc = JKmMkk{l)xk + Rf + Nfc (3) 



K 

,{SD) 



Rl - E ^P^Pk^^k^{l)Xi (4) 

In these equations, Y, H, R and N are all M x 1 vectors. R denotes the inter-pair interference. 
Xk is the transmitted symbols for the fcth pair, which has zero mean and unit variance. Pk is the 
allocated power for Xk- is a white Gaussian noise vector with covariance (j^Ia/. 

2) 2 X M STBC: In the 2 x M Space-Time Block Coding (STBC) scenario we use an Alamouti 
Code [fTSl with two transmit antennas. We use Xk,m,n to denote the transmitted symbol at mth 
antenna and nth time slot of the kth pair. Using the Alamouti Code, we have Xk,i,2 = —xl 2 1 
and Xk,2,2 = xlii, {1 < m < 2,1 < n < 2). Here Xk,i,i and Xk,2,i are independent symbols 
with zero mean and unit variance, and (■)* is the complex conjugation. Pa,(1) and Pa;(2) are 
allocated power for Xk,i,i and Xk,2,i, respectively. Assume the received vectors corresponding to 
these two time slots are Yk^i and 2, as well as the noise vectors Na:,i and Nfc,2- Define new 
vectors by = [Yk,i]"^l2,], ^ki,i = [Hfei(l); Hfci(2)*; ], Hki,2 = [Hfci(2); -Hh(1)*; ] and 
Nfc = [Nfc 1; 2! ]• Consequently, the received signal at Rx Node k can be denoted as: 



Pk{l)pkkiikk,iXk,i,i + \/Pki2)pkkiikk,2Xk,2,i + RI'^^'^'^^ + Nfc (5) 

K 2 

, (STBC) 



Rl ^= E Ev^^(0pfc3fc,,„x,,z,i (6) 
i=i,i^k 1=1 

3) M X M Transmit Beamforming: The Beamforming method outlined in [fT6l. [[T7| is also 
employed in this paper. M transmit antennas are adopted, and an M x 1 weight vector is 
applied to these transmit antennas. Let Hkk = [Hfcfc(l), Hfcfc(2), . . . , HkkiM)], then according to 
[fT6l . [fTTl . Ufc is the eigenvector corresponding to the largest eigenvalue of the matrix H^^H^fc, 
where (■)^ denotes conjugate transpose. We assume that uj^Uk = 1 and Xk has zero mean and 
unit variance. The transmit power for Xk is Pk. Then the received signal under this method is 
represented as: 

Yfc = \l Pkpkk^kkVLkXk + R^^*"""'^ + Nfc (7) 

K 

I (Beam.) 



Rfc = E \l PiPkitikiViiXi (8) 

i=l,i^k 
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C. MMSE Solution 

For all three scenarios listed in the previous section we will use the Minimum Mean Squared 
Error (MMSE) solution at the receiver |[T8l, [[l9l, lED to arrive at the receive MIMO antenna 
weights. Here we only present the MMSE solution for the 1 x M Receive Diversity case. 
However, the MMSE solution can also be derived for the 2 x M STBC and M x M Transmit 
Beamforming scenarios by slight modification (the associated results are given in Appendix I). 

For the 1 x M Receive Diversity scenario, the estimate of the transmitted stream at the 
receiver is given by Q where the vector is the MMSE solution. 



Xk = wf Yfc (9) 



The corresponding SINR at the receiver is: 



(10) 



^iSD){k) = PiPfc^Hfci(l)H^(l) + a^lM (11) 

l<i<K,i^k 

With the constraint that wj^Hk^i^) = 1, optimal linear vector w^. in the MMSE solution is 
given in El, EH, [EH 

^ $^^)(fc)H,,(l) 

The SINR with MMSE solution is: 

ffc = PfcPfcfeHf,(l)<l>(-^)(A;)H,fc(l) (13) 

Finally, we assume a QoS requirement at the receiver. Consider an SINR threshold 7^, for 
the kth transceiver pair, they can be correctly received if and only if the received SINR is not 
lower than 77'. The threshold 72- represents the QoS requirement and is being used here in the 
same manner as the Physical Model in [|20ll . 
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III. The Definition of Concurrent Links 

We assume that all Tx and Rx Nodes are mobile, and their locations are randomly distributed 
and varying with time. In the following discussion, a specific scenario means one specific 
realization of node locations in the region, and the channel responses between them. For a specific 
scenario, given the allocated power in each pair, the SINR for each pair can be evaluated using 
the results of Section n. Based on these SINR results, we first provide the following definition. 

Definition 1 (Link): Consider a specific scenario with K transceiver pairs. The transmit power 
for every pair is constrained by Pt. The fcth pair (1 < A; < i^) is called a link if and only if it 
satisfies the following conditions: 

For 1 X M Receive Diversity and M x M Transmit Beamforming: 

ffe > 7T and Pk < Pt (14) 

For Space Time Block Coding: 

tkil) > 7t, 1 < i < 2 and P^il) < Pt (15) 

1=1 

Assume a specific scenario with K transceiver pairs, labeled from 1 to i^' and denoted by a 
setC/ = {1,2, A pair set C/p C [/ is a feasible pair set if there exists a power allocation 
for all transmitters in Up such that Definition 1 holds for all transceiver pairs in Up (i.e., the 
QoS constraint is satisfied at all receivers subject to the maximum Pt constraint). Then all the 
pairs in C/p are named Concurrent Links. 

Let us denote the number of pairs in Up as: \Up\. Then the MNCL is calculated by the 
following optimization problem: 



max \Up\ 
s.t. Up CU and Up is feasible (16) 

For notational simplicity, the result of the above optimization problem is denoted as A^max(^, M). 
Note that several different feasible sets Up may produce the same Nais,^{K, M). We then define 
Average MNCL as the expectation of iVmax(^, M), averaged over random locations and random 
channel responses. The Average MNCL will be denoted as C{K, M) in this paper. 
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The optimization in equation (fT6l ) can be divided into two steps: 1) to examine whether a 
pair vector Up is feasible; and 2) to search for the feasible pair set U p with the MNCL. In this 
paper, the first step is solved by an optimal power allocation process presented in Section IV, 
while the second step, referred to as optimal link selection, is solved in Section V. 

IV. Optimal Power Allocation 

Optimal power allocation is an important factor when deriving the bound on the MNCL. In our 
formulation it is used to decide whether a given pair set Up is feasible per Definition 1. In this 
section, the algorithm for optimal power allocation is presented from conventional power control 
techniques ETTl - lllBl . Without loss of generality, only the 1 x M Receive Diversity scenario is 
considered in this section. However, the proposed algorithm can also be applied to the 2 x M 
STBC and M x M Transmit Beamforming scenarios after slight modification. 

A. Iterative Power Control 

Consider K transceiver pairs in the network using 1 x M Receive Diversity. The allocated 
power for each pair is stacked into a vector P, which is defined as the power vector and given 
by: 

P = [Fi,P2,...,Pi^] (17) 
Since we only focus on pairs in Up, we have the following constraint: 

for ^ f/p, Pfc = 0, l<k<K (18) 

Let Wfc be the MMSE solution in equation (fT2l) . We define the following iteration equation: 

For keUp 

p,,. _ ^^ Cr' {w.,P"} + ai,wfw> ^^^^ 

Pkk 

cf^^{w,,P'^}= Ofc.G'{w,,HH(l)} (20) 

G {wfc, Hfc,(l)} = wf Hfc,(l)H^(l)wfc (21) 
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For l<k<K andk ^ Up 

P^+' = (22) 
where n denotes the nth iteration. For simplicity, the above iteration is denoted as 

P"+i = m (P") (23) 

Here P" is the power vector for the nth iteration. 

We define the fixed point of mapping as the power vector P satisfying P = m(P). The 
following theorem holds for the iterative equation (|23l) . which is used to verify the existence of 
optimal power allocation in this paper. 

Theorem 1: Given K transceiver pairs and a specific pair set Up. If Up is feasible, then a 
unique fixed point of mapping, P, exists that satisfies P = m(P) and Pk < Pr.^k E Up. 
Furthermore, corresponding to the unique power vector P is a unique receive weight vector 
given by the MMSE solution. 

Proof: Please refer to Appendix II. Note that the converse proposition of this theorem holds 
obviously. That is, if there exists P = m(P) and i\ < -Pt,VA; E Up, then Up is feasible. ■ 

B. Decision Criteria 

Theorem 1 suggests that a feasible pair set can be identified if and only if a fixed point of 
mapping for the transmit power levels exists and the power constraints are met. Moreover, f2l\ 
suggests that the power control algorithm in (|23|) is guaranteed to converge to a fixed point of 
mapping even when starting from an arbitrary initial power vector P°. This suggests a rather 
straight forward approach for determining the feasibility of a given pair set Up as captured by 
the following three criteria. 

Criterion 1: Assume that the iteration process starts from initial condition P° = 0. In each 
iteration, if the power in any transceiver pair exceeds the power constraint Pt, then Up is not 
feasible. 

The proof of Criterion 1 is referred to Theorem 2 in [21]. 

Criterion 2: Assume that the iteration process starts from arbitrary condition P°. In each 
iteration, if P" is feasible, then these pairs can be supported simultaneously. 
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Here, saying P" is feasible means that by using P", for 1 < k < K, we can have > 
and PI' < Pt- This criterion is straightforward and it is helpful in reducing the iterations in the 

decision process. 

Criterion 3: After /max number of iterations, if the conditions in Criterion 2 have not yet 
been satisfied for pairs in Up, then we declare that Up is not feasible. 

Criterion 3 guarantees that the decision can be made within a finite iteration number /max- 
Using this criterion, some feasible pair sets may be missed. However, proper selection of /max 
can drive the probability of missing a feasible pair set to be arbitrarily close to zero. 

Using the above criteria we now outline a sequential procedure by which we can automatically 
determine if a set ?7p is feasible. The steps are as follows: 

Iterative Determination of Feasibility (IDF) 

1) Given pair set Up. Initialize n — Q and P° = 0; 

2) Iterate by P'^+i — m (P"). Using the updated power vector, P"+^, calculate the SINR for 
each pair under the MMSE solution; 

3) liUp is not feasible by Criterion 1, then go to step 6; else, go to step 4; 

4) If Up is feasible by Criterion 2, then go to step 6; else, go to step 5; 

5) If Up is not feasible by Criterion 3, then go to step 6; else, n = n + 1, go to step 2. 

6) The feasibility of Up is determined, stop. 

V. Finding Up with Maximum Number of Concurrent Links 

In the previous section, we derived the method to determine if a pair set is feasible. In this 
section we describe how to find the feasible pair set with the MNCL. Naturally, a brute force 
search can be employed for this problem. However, in order to reduce the search space, we first 
characterize the property of the feasible pair set. 

Consider a total of K transceiver pairs, labeled 1 to K. Define U = {1,2,..., K} and Up is 
a given pair set. 

Theorem 2: For two pair sets Up^2 Q Up^i, if Up^i is feasible, then C/p2 is also feasible. 
Proof: Consider the 1 x M Receive Diversity scenario. Let (?7p,i — [^p,2) denote the pairs 
in f/p i but not in f/p,2- If f^p,i is feasible with the associated power vector P, then for the pairs 
in Up,2, keep the power P^ and linear vector unchanged. Meanwhile, shut down the pairs in 
{Up,i — Up^2)- Here, shutting down means setting the corresponding transmit power to zero. As 
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a result, there are less interference in ^(^sD)ik) of equation (10), and conditions in Definition 1 
are satisfied for the pairs in f/p2- Thus f/p2 is feasible and the above theorem holds. The proof 
is extendable for STBC and Beamforming methods as well. ■ 

Theorem 2 shows that the backtracking formulation in |[8l, [|24ll can be adopted to solve 
the problem of finding the Up with the MNCL. We start this formulation with an empty pair 
set. At every level, each feasible subset is expanded by including one more pair, constructing 
new subsets to be validated by the IDF iterations presented in Section IV (namely, forward 
search). If one subset becomes unfeasible, the algorithm backtracks by removing the trailing 
pair from the subset, and then proceeds by expanding the subset with alternative pairs (namely, 
backward search). Specifically, we use a depth-first search strategy to execute this procedure. 
An example with 4 pairs are illustrated in Fig. 2, where all feasible pair sets are outlined. 
Here the initial feasible subsets are {{1}, {2}, {3}, {4}}, while the final feasible subsets are 
{{1, 2, 3}, {1, 3}, {2, 3}, {2, 4}, {3, 4}, {4}}, and the MNCL is 3. 

Let Up^search dcuotc the pair set candidate during the search, and |f/p,search| is the number 
of elements in Up^search- Meanwhile, use direction = 1 to represent the forward search, and 
direction = the backward search. The pseudocode of the proposed algorithm is described in 
the following. 

Backtracking-based Optimal Link Selection (BOLS) 

1) Given K pairs (labeled from 1 to K) in the network and M receive antennas per pair. 
Initialize by Up^search = {1}, N„,^^{K,M) = 0. 

2) If \Up^search\ > N^^^{K,M), thcu: 

a) Use the IDF process in Section IV to examine whether U p^search is feasible. If U p^search 
is feasible, then set NmaxiK,M) = \Up^search\ and direction = 1; otherwise, set 
direction = 0. 

3) If \Up^search\ < N^!^^{K,M), thcu sct dircctiou = 1. 

4) Update the pair set candidate Up^search- 

a) If Up^search = {K} , thcu all the search have been done, return A^max(-^5 and stop. 

b) If Up^search 7^ {K} , thcu Up^search = PairSet_Gen [Up^search, direction] and go to 
step 2. 

In step 3, we only examine the pair set Up^search in which \Up^search\ is larger than the current 
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value of N^s.y.{K,M). In other words, if \Up^search\ < Njjig,^{K, M), instead of employing the 
IDF procedure, we set direction = 1 and go to step 4 directly. Next, the function PairSet_Gen 
is described as follows: 
Function PairSet_Gen 

1) Input Up^search and direction. Find fcmax = max k, which represnts the pair with the 

^^^P, search 

maximum index in Up^search- 

2) If A;max < K, then: 

a) If direction = 1, add new element in Up^search by Up^search = Up^search U {fcmax + !}• 
Return Up^search and stop. 

b) If direction = 0, update the element fcmax in Up^search by fcmax = ^max + 1- Return 

Up^search and StOp. 

3) If A^max = K and \Up^search\ > 1, then; 

a) Delete /C^ax from Up^search by f/piearc/i = Up,search-{kma^}- Find fcmax = HiaX k, 

search 

update the element Ai^ax in f/pfiearch by fc^ax = A^max + 1- 

b) Set Up 

search = ^P^searc/i- RctUm Up^search and StOp. 

4) If /cmax = K and |f/p,searc/i| = 1, then all the search have been done. Return Up^search and 
stop. 

Finally, averaging A^max(-^, M) among random locations and random channel responses, we 
obtain the C{K, M). 

VI. Simulation Results 

The simulation setup randomly distributes the nodes, in accordance to a uniform distribution, 
within a disk of radius 100 meters. Numerical results are averaged over 1000 Monte Carlo 
simulations of independent realizations of node topology and channel response. We assume 
QPSK modulation in the baseband, and the desired SINR threshold, 72-, is lOdB. This corresponds 
to an uncoded BER of less than le-3 (see Table 6.1 in lfT4l '). The remaining parameters are listed 
in Table I. 

A. Comparison of MIMO Techniques 

In this section we evaluate the performance of the proposed algorithms, and use the MNCL as a 
metric to compare different MIMO methods. The three MIMO methods outlined in Section II are 
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compared, namely, 1 x M Receive Diversity, 2 x M STBC, and M x M Transmit Beamforming, 
all with MMSE decoding. 

We assume up to 15 transceiver pairs, and vary the number of receive antennas from 1 to 4. 
Results for the three different MIMO methods are reported in Fig. 3 and Fig. 4. We observe that 
Transmit Beamforming has the best performance due to the fact that it explores Channel State 
Information (CSI) at the transmitter. On the other hand, STBC has a worse performance than 
Receive Diversity. Note that STBC is conventionally designed to cope with the white Gaussian 
noise [flSl . but in this study we use MIMO system to combat the colored interference signals 
(equation (fTT|) ). and our results indicate that Receive Diversity outperforms STBC in this case. 
Finally, performance of conventional SISO system corresponds to that of Receive Diversity with 
M = 1. We observe that compared to a SISO system, as much as 3x improvement in MNCL can 
be achieved by using M = 4 receive antennas (compare the highlighted values for K = 5, 10, 15 
in Fig. 3 with those in Fig. 4). 

Next, we study the convergence of the estimate for the MNCL. We assume a total of 12 
transceiver pairs in the network, with the number of receive antennas fixed 4. As previously 
mentioned, we create a realization by distributing the pairs in a 100m radius disk and randomly 
generating the channel responses among them. In Fig. 5, we plot the average number of MNCL 
versus the number of realizations. We find that the simulation results become convergent after 
around 500 realizations. Note that we used 1000 independent realizations in these simulations, 
which is large enough to yield a precise estimate for the performance of the proposed algorithm. 

Finally, we compare the proposed backtracking strategy (Section V) with the brute force 
search. We have verified that, for I < K < 15, the proposed backtracking method has exactly 
the same MNCL as the brute force search (the results are omitted due to space limit). Here we 
focus on the search complexity of these two search methods, that is, the number of times the 
IDF iterations were called. The average result over 1000 Monte Carlo simulations is shown in 
Fig. 6. Note that the associated complexity for the brute force method is 2^, where K is the 
number of pairs in the network. Compared with the brute force search, we observe a significant 
reduction in complexity when the backtracking scheme is used. This verifies the efficacy of our 
proposed scheme. 
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B. Impact of Simulation Parameters 

Using the formulation developed in this work, we can investigate the impact of different 
parameters, such as the number of transceiver pairs, the number of receive antennas, and the 
maximum allowable transmit power on the MNCLs. We first look at the number of pairs in the 
network. Results in Fig. 3 and Fig. 4 have already shown that the MNCL increases substantially 
with the number of pairs. Actually, we can use a two-stage linear equation to approximate the 
result of the MNCL. The equation is presented in (|24|) . and the approximation result for Transmit 
Beamforming is demonstrated in Fig. 7 and Table II. 

For I < K <lh (ai, 6i, 02, 62 are parameters to be fitted): 



In the above equation, when 1 < K < M, the improvement is mainly from the diversity gain; 
while when M + 1 < K < 15, the improvement is from multi-user diversity (as the number of 
pairs wanting to communicate increases the likelihood of choosing a subset of these pairs that 
exhibit good interference properties increases). The results in Fig. 3 and Fig. 4 show that the 
MNCL is larger than M when the number of pairs K is increased beyond M. This observation 
highlights some extra gains in MNCL to be exploited by the MAC protocol. 

Next, we focus on the impact of the number of receive antennas. We assume 12 transceiver 
pairs, and the number of receive antennas M varies from 1 to 8. Simulation results are shown 
in Fig. 8. We see that in this simulation, the MNCL is improved dramatically with the number 
of receive antennas. Again, we find that with 4 receive antennas, the improvement is around 3 
times relative to the SISO system. Furthermore, we observe that saturation starts to set in when 
the number of receive antennas is greater than 5. With 8 receive antennas the gain relative to a 
SISO system is 4x. 

Lastly, we explore the impact of the maximum allowable transmit power Pt. Naturally, higher 
transmit power is useful to combat path loss and increase transmission range. However, here 
we analyze the impact of under the assumption that all the nodes are distributed in a fixed 
disk with radius 100m. This assumption is reasonable which corresponds to a geographically 
constrained area with many potential transmission pairs (e.g., classroom, conference room). We 
assume a total of 10 pairs in the network, with 4 receive antennas per pair, and consider different 
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power constraints from -20dBm to 50dBm. The results are depicted in Fig. 9. We observe that 
for each MIMO method, initially, the MNCL increases with higher power constraint. However, 

when Pt is beyond lOdBm, the improvement becomes diminishing. The reason is that co-channel 
interference among pairs has dominated the network, and increasing Pt can not mitigate these 
interfering signals. 

VII. CONCLUSION 

In this paper, we have identified the bound on the maximum number of concurrent links 
(MNCL) in MIMO Ad Hoc networks subject to a minimum QoS (SINR in our case) for each 
link. This number is derived by considering the practical factors of wireless channel, transmit 
power allocation, link selection and MIMO transceiver algorithms in a realistic MIMO system. 
We employ an iterative algorithm to examine the existence of an optimal power allocation that 
guarantees QoS for all pairs in a given set. Based on this iterative algorithm, we proposed a 
backtracking search algorithm to select the optimal subset of pairs, constituting the MNCL. 

Extensive simulations were conducted to verify the efficacy of the proposed algorithms and 
evaluate the impact from different parameters. The results show a 3x improvement in MNCL 
with a 4 element antenna system relative to a SISO system. For the parameters simulated, 
diminishing improvement is observed when the number of antennas is increased beyond 5, and 
when the maximum transmit power is increased beyond lOdBm. 

Appendix A 

In this appendix, we show the MMSE solution for STBC and Transmit Beamforming methods. 
For STBC, define that the MMSE solutions for Xk,i,i and 0:^^2,1 are w^^i and 2, respectively. 
Then we can have: 

^{STBC,l)i^)'^kk,l 
Wfe,l = ' — — = — (25) 

^kk,l^(STBC,l) (^)Hfefc,l 
^(SrBC,2)(^)Hfefc,2 

Wfe,2 = _H — — = — (26) 

Hfcfe,2^(5TBC,2)(^)Hfcfe,2 

K 2 ^ ^ 

^{STBC,i){k) = Z]^i(OPfciHfci,«Hfci; + Pfc(2)pfcfcHfefc,2H;i,fc,2 + c^atIm (27) 

i=l,i^k 1=1 
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K 2 

^{STBC,2)ik) = J2 Pii^)PkMki,i^i,i + Pk{l)PkkH.kkA^fk,i + ^a'Im (28) 

i=l,i^k 1=1 

For Transmit Beamforming, the MMSE solution is: 

K 

^(BF){k)= i^iPfci(HfciUfc)(Hfe,Ufc)^ + a^lM (30) 

Appendix B 
Proof of Theorem 1 

Note that we have assumed in (fTSi) that {Pk = 0,k ^ Up}. The iteration equation (|23l) can 
be proved to be a standard function ll23l for {Pk,k E Up} via similar manner in [[2T|. Now 
assume that with power vector P, Up is feasible. Set the initial power vector P° = P, and run 
the iteration process by P"+i = m(P"). According to Lemma 1 in ll23l . P" will convergence 
to the fixed point of mapping, which is P = m(P). With Mono tonicity property in standard 
function, for 1 < k < K, we have Pk < Pj^ < Pt- Finally, it has been proved in ETI that 
power and weight vectors corresponding to the fixed point of mapping are all unique. Thus the 
theorem follows. 
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TABLE I 
Parameters in the Simulations 



Parameter 


Value 


Path Loss 


do ~ Im 




Lp(do) = 46dB 




Exponent factor a = 3 


Noise Power Spectral Density 


rir, = -174dBm/Hz 


Noise Figure 


F„ = 4dB 


Bandwidth 


11= IMHz 


Maximum Transmit Power 


Pt = 20dBm 


SINR threshold 


7T = lOdB 



TABLE n 

Parameters in the Approximation Results of Average MNCL with Transmit Beamforming 



M 


1 


2 


3 


4 


ai 





0.0040 


-0.0003 


0.0030 


bi 


0.9680 


0.9950 


0.9990 


0.9956 


£12 


0.9961 


1.9587 


2.6045 


3.5086 


b2 


0.1084 


0.2121 


0.3013 


0.3171 



Rx Node O Tx Node 

2 / \ 2 Desired 

Communication 



— 

Node A ^\ / ►Q Rx Nc 

1 y y\ A 1 



Tx Node i^^Z Rx Node 

1 y ^\ ' Interference 

Rx Node 1 / / Tx Node ► 

3 3 



Fig. 1. Illustration of the considered network. SoUd lines are the desired communication, while dashed Unes are the interference. 
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Fig. 2. Tree structure for backtracking search. 
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Fig. 3. Average MNCL under different MIMO teciiniques. M is from 1 to 3. 
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Fig. 4. Average MNCL under different MIMO techniques. M is set as 4. 
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Fig. 5. Estimate results for the MNCL versus the number of independent realizations. Assume total 12 transceiver pairs and 
M = 4. 
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Fig. 6. Average number of times the IDF procedure was called versus the number of transceiver pairs for M — 4. 
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Fig. 7. Simulation and approximation results for average MNCL with Transmit Beamforming. Assume 12 pairs equipped with 
4 receive antennas. 
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Fig. 8. Results under different numbers of receive antennas. Assume total 12 pairs. Up to 8 receive antennas are simulated. 




Fig. 9. Average MNCL versus different maximum allowable transmit power per pair. Assume a total of 10 pairs, Af is fixed 
as 4. 
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